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The phase diagram of the one-dimensional extended Hubbard model at half-filling is investigated 
by a weak coupling renormalization group method applicable beyond the usual continuum limit for 
the electron spectrum and coupling constants. We analyze the influence of irrelevant momentum 
dependent interactions on asymptotic properties of the correlation functions and the nature of 
dominant phases for the lattice model under study. 



PACS numbers: 



I. INTRODUCTION 



The application of the renormalization group (RG) and 
boson representation methods to one-dimensional (ID) 
models of interacting electrons have provided over the 
last four decades considerable insight into the nature 
of correlations in low dimensional systems 1 -^. This has 
been largely achieved by treating the models in the con- 
tinuum field-theory limit, corresponding to the so-called 
weak coupling ID electron gas (EG) model. There are 
notable exceptions, however, where the 1D-EG model 
clearly fails to reveal the nature of correlations at long 
distance. These situations are likely to occur in models 
for which lattice effects, albeit related to irrelevant terms 
in the RG sense, do affect the asymptotic behavior of 
electronic correlations and then the nature of the ground 
state. 

A well documented case is encountered in the extended 
Hubbard model at half-filling, which is defined on a lat- 
tice in terms of intersite hopping and the on-site and 
nearest-neighbor sites couplings U and V. On numeri- 
cal side, exact diagonalization^, quantum Monte Carlo^ 
and density-matrix renormalization group analysis^ have 
established the incursion of a bond order-wave (BOW) 
state over a finite region of the phase diagram surround- 
ing the line U = 2V > 0, a result at variance with the 
spin density-wave (SDW) to charge density-wave (CDW) 
transition found in the theory of the ID EG 1 ^^. 

Using perturbation theory arguments, Tsuchiizu and 
Furasakiii showed how high-energy or short-distance de- 
grees of freedom can modify the initial conditions of an 
effective low energy continuum theory and favor the oc- 
currence of a BOW phase that enfolds the U — 2V > 
line in weak coupling. The influence of the lattice on the 
nature of the ordered phase in this region of the phase 
diagram has been investigated by Tarn et al.jH using the 
functional RG method. The scaling transformation of 
interactions, which in this framework gather both their 
marginal and momentum dependent parts, was obtained 
for a tight-binding electron spectrum in a finite momen- 
tum space. The CDW/SDW degeneracy that takes place 
at U = 2V in the continuum limit, is thus lifted and a 
BOW state stabilized over a portion of the phase dia- 



gram that grows in size with increasing U , consistently 
with numerical calculations at weak coupling. The func- 
tional RG method, however, tells us not as much about 
the structure of leading irrelevant terms and how these 
modify scaling and the nature of ordered states of the 
continuum theory. 

In this work we address this issue from a different 
perspective that generalizes the weak coupling momen- 
tum shell Kadanoff- Wilson (K-W) RG method to lattice 
model s 13 ' 14 . The proposed approach exceeds the limi- 
tations of the continuum approximation and takes into 
account the tight-binding structure of the spectrum and 
its impact on the scaling transformation of both local and 
momentum dependent interactions of the extended Hub- 
bard models. The latter couplings, though irrelevant, 
are found to affect the flow of the former interactions. 
A modification of certain portions of the phase diagram 
follows; in particular, the BOW phase is found to insert 
in a finite region near the U = 2V > line, in agreement 
with the results of numerical calculations. 

In Sec. II, we introduce the model and set out the basic 
steps of the momentum shell RG transformation for the 
partition function. In Sec. Ill, the RG flow equations for 
the coupling constants and the most singular response 
functions are analyzed at the one-loop level and different 
U and V. The phase diagram is mapped out in weak 
coupling. We conclude in Sec. IV. 



II. THE EXTENDED HUBBARD MODEL AND 
THE RENORMALIZATION GROUP 
FORMULATION 



A. The model 

We consider the extended Hubbard Hamiltonian for a 
one-dimensional lattice, 



(1) 
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where t is is the hopping integral, n ijCr = cJ CT c iiCr is the ,\., and = n^+ni^. In Fourier space the Hamiltonian 
occupation number on site i for the spin orientation a =\ can be written in the form 



H = tkclk^c p .k.* +- ^ (.9i + 25ism 2 -jc^ fci+9+2fcF ^c^ fe2 _ g _ 2feF CT2 c +!fe2i(T2 c_ !fcl ^ 1 

p.k.cr {k,q,a} 

+ ^ E (.92 + 2.9 2 sm 2 |)c 1 f!fci+(J(Ti cL ifc2 _ 9(T2 c_, fe2 , (T2 c + ^ l , (Tl 

{k,q,o} 

+ JZ (.93 + 2.93 sin 2 |) (c^ ifcl+9+afcF) ^c^ fe _ ff _ 2fcji . +0)OJ1 c_,* a ,^c_,fc 1 , ffl +H.c.) ( 2 ) 

{fc,g,o-} 

+ ^ Z (.94 + 2.94 sin 2 |) c^ M+q ^J +M _ q ^ 2 C +M ^c +M , ai 

{k,q,a} 

+ JZ Z (-94 + 254sin 2 |) cL fci+g (Ti cL !fc2 _ gi(T2 c_, fc2 , . 2 c_ :fcli<Tl , 

{fe, <J,er} 



where = — 2tcosk is the tight-binding spectrum, and 
v = 2t is the bare Fermi velocity; the Fermi points are 
kp = ±-| at half- filling (here the lattice constant has 
been set to unity, and h = 1 = feg). By analogy with the 
'g-ology' description of interactions, we have proceeded 
to the splitting of the U and V interaction terms into 
couplings for right (p — +, k > 0) and left (p = — , 
k < 0) moving electrons. We thus obtain momentum in- 
dependent (local) as well as momentum dependent (non 
local) couplings, denoted in @ by Qi=\...\ and g\, ..4, 
respectively. The pairs of couplings for backscattering 
(fl'iiSi) and Umklapp (33,33) have the bare amplitudes 
.91,3 = U — 2V, 31,3 = 2V, whereas 32.4 = U + 2V and 
.92,4 = —2V stand for the amplitudes for the forward 
scattering between opposite (32,32) and parallel (34,34) 
k electrons. 

The information about the lattice in @ is present by 
the use of the tight binding spectrum e k for k e [— n, ir] 
in the Brillouin zone and in the momentum dependent 
couplings 3^ sin 2 9/2. In the continuum limit, the latter 
amplitudes vanish when evaluated at zero momentum 
transfer, while the spectrum e k — > e p (k) ~ v(pk — kp) 
is taken as linear around each Fermi points. One thus 
recovers the standard electron gas formulation of the ex- 
tended Hubbard model*—. 



B. The renormalization group transformation 

We write the partition function Z = Tr e~^ H as a func- 
tional integral 

Z = J J Dip*T)ip e s W (3) 

over anticommuting Grassmann fields ip(*'. The action 
S[ip*,ip] = So[ip*,tp] + Si[ip*,tp] consists of a free and an 
interacting parts. In the Fourier-Matsubara space, the 
former part So[ip*,tp] reads 

Soir,^} = E K^T^lA^vAk), (4) 

p,k,cr 

where 

G° p (k) = [iuj n -e k ]-\ (5) 

is the free electron propagator. Here k = {k,uj n ) and 
Lu n = (2n + 1)tvT is the fermion Mastubara frequency. 
The interacting part is given by 
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Si[^*,tp] = --=■ X { 9l + 2 5i sin2 |) V'+, ( T 1 ( fc i +90 + 9)^1,^(^2 -qo~ 3)^+,<t 2 (h)^-,^ (fci) 



{k,q,a} 



T 



,2 9 



1 zZ \92 + 2g2^ 2 ±)r+^ 1 {h+qW-,a 2 (k2-q)^^ 



{k,q,a} 



T 

2L 

T 
2L 



X (-93 + 2 53 sin 2 |) (ip+ t<T1 (fci +q + g#+,<r a ( fc 2 - 9o ~ 9 + G)i>-,* a (fo)V>-<n (*i) + c.c.) (g) 

{fc,g,o-} 

X] (-94 + 2.g 4 sin 2 |) (*i + e)^,^ - g)^+,a 2 (fe)^,^ (fci) 

{fc,g,o-} 



57 zC (f4 + 2^4 sin 2 I) + <l)ip-,a 3 (k2 - q)^-,a 2 (k 2 )^-,a 1 (ki), 



{k,q,a} 



r 



where q = (g,w m ), w m = 27rmT, go = (2fcF,0); here 
G = (4fcj?,0) is a reciprocal lattice vector that enters in 
the definition of Umklapp scattering at half-filling. 

The momentum shell K-W RG transformation is based 
upon the recursive application of the two following steps 
for the partition function. In the first step, a partial 
trace of Z over outer shell electronic degrees of freedom 
denoted by ipp >t r(k,u} n ), is carried out at all u> n and spin 
cr. The outer momentum shell is denned by the intervals 
of momentum 



k € [0, hp — ko/s[ U + fco/s, 7r], p = + 
S ]— hp + ko/s, 0] U [—w,—kp — ko/s[, p = 



(7) 



above and below the Fermi level for each branch p. Here 
fco = 7r/2 is a cutoff wave vector denned with respect to 
the Fermi points (±fep ± fco = ±tt), and s = e > 1 is 
the momentum scaling factor for < 1. The second 
step consists in the rescaling of the momentum distance 
from the Fermi points we call Sk; this gives k! = ±k F + 
sSk, which restores the initial cutoff fco = tt/2 of the 
lattice model. 

The two recursive steps of the RG transformation can 
be expressed as 



Z 



x 



x e £?=i -Si,* [*/>*, </>,'</>*,'/'] 

£)lp*Dlp e S '[ , / , *^]i + <S'7,2>0, c + I<(S'/,2) 2 >o,c+-' 



(8) 



where Sj , is the interacting part of the action with i = 
1, . . . , 4, ip fields in the outer momentum shell. The outer 
shell statistical averages (■■■■)o, c over the variables 
are performed with respect to So [V>* , V»] . These averages 
correspond to the sum of all connected diagrams with 



even number of external fields ip*,ip pertaining to the 
inner momentum shell (<) degrees of freedom, which are 
kept fixed in the partial trace operation. At the one-loop 
level the partial trace and rescaling lead to the recursion 
relations 



S [ip*,ip]e +d e = S [%jj*,ip]e + (Si, 2 )o, c + 



(9) 



ifi—tCfifi' 



(10) 



for the free and interacting parts of the action. 

Following the momentum rescaling, the inner shell 
fields if) are rescaled by the factor Cs^ 2 , which can be 
derived from a dimensional analysis of the parameters 
that define the bare action So- Thus assuming that the 
rescaling of the tight-binding spectrum is of the form 
e' k , = Cs^fc, by taking k' — ±kp + sSk, one gets in the 
limit di -> 



Sk cot Sk 



(11) 



which can be expressed in the form s v compatible with 
an iterative transformation in renormalization group. At 
variance with the usual case, however, the scaling di- 
mension y is here fc dependent. Thus at either the edge 
or the bottom of the band where the group velocity 
vanishes, Cs(±fco) s° and is dimensionless. It is 
only when the Fermi points is approached in the limit 
Sk —> 0, that C s (Sk) — >■ s 1 and the result of the con- 
tinuum limit for a linear spectrum is recovered^. This 
also indicates that repetition of rescaling turns down the 
curvature of the band, which continuously evolves to- 
ward a linear shape. Since uj n or the temperature T 
enters on the same footing as efe in the inverse propaga- 
tor [G ] -1 , the temperature then transforms according 
to T' = £ s (Sk)T. Now referring to the form of Sq in 
0, this yields the transformation assumed above for the 
field, namely = C7 1/2 {Sk)ipW. 
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When applied to the interacting part Si of the the 
action, the above relations for the field and temperature, 
combined to the shrinking of the number sites V = L/s 
under rescaling, will impose the following fc-dependcnt 
transformations of interactions 



0{g 2 )) t 
■0(gg))t 



-l-\-Sk cot 5k 



-1 — Sk esc Sk 



(12) 
(13) 

It ensues that for Sk — > ±fco, we have g[ — > 5iS _1 , and 
the local couplings are then irrelevant instead of being 
marginal variables near the bottom or the edge of the 
band. At the approach of the Fermi level, when Sk — > 0, 
we have g[ = gi and the dimensionless or marginal char- 
acter of the local interactions of the electron gas model 
is retrieved^^. In the same way, the non local terms 
transform according to g[ — giS" 1- ™/ 2 at the boundaries 
or the bottom of the band and are therefore strongly ir- 
relevant. In the limit Sk — > near the Fermi points, 
9i — gi s 2 i which corresponds to the usual negative bare 
scaling dimension of nearest-neighbor couplings of the 
continuum theory^^. 



C. The Fermi velocity and coupling constant flow 
equations 

We now proceed to the partial trace operation that 
defines the first step of the renormalization group trans- 
formation ©. At the one- loop level, this amounts to 
evaluate the outer shell statistical averages (<S7,2)o c an d 
(Sj 2 )q c of the recursion relations © and (fTU)) . The for- 
mer contribution (S7,2)o.c i s composed of Hartree and 
Fock self-energy corrections. In these, enter k indepen- 
dent or constant terms that correct the chemical poten- 
tial, a quantity that can be simply redefined to keep the 
filling of the band constant. These terms can be safely 
ignored. The presence of non local interactions give rise 
to momentum dependent Fock terms, which at the step 
£ of the iterative RG procedure read 



( s 1,2)0,, 



L{£) 



££[^)G%(F)-s 4 (l)G£(fc')] 



k k> 



x cos(kp + 5k)ip*(k)ipp(k), 



(14) 

where the slashed summation contains an integration 
over kl in the outer momentum shell interval ([7|1 at a 
given p. The Fock terms contribute to the renormaliza- 
tion of the spectrum, that is the Fermi velocity. Carrying 
the k' summation, one gets the flow equation for the ve- 
locity, 

d e lnv(£) = |(<? 4 W -9i(t)) tanh[w(f)sin(5fc £ /2T], (15) 

where Ski = koe~ e and the couplings g = g/wv(£) are 
henceforth taken as normalized by the scale dependent 
Fermi velocity v{£). 



The recursion relations (|12l) for the local normalized 
couplings <?(= g/irv{£)) are obtained from the outer 
shell contractions {Sj 2 )o jC in the logarithmically singu- 
lar Cooper (electron-electron) and Peierls (2kp electron- 
hole) channels. Their insertion in (fT2"]) . leads after rescal- 
ing to the recursion relations 

9i = [.9i + (Si + hh)!p + .9i (.92 + hVc }s~ fa (16) 
.9 2 = [h+ {h+h) 2 ic + {gz+hfip\s' h (17) 



93 = [h + (<72 + 92X2.93 + h)!p 



{91 +9i)(d3 -g3)ip]s 



-f g 



94 s 



(18) 
(19) 



where f g = 1 — Ski cot Skg + di In v(£) , which contains the 
rescaling exponent of (TT2"j) . and the correction due to the 
normalization from the scale dependent Fermi velocity. 
We note that the one-loop level, there is no logarithmic 
correction to the forward scattering amplitude 54. The 
outer shell Cooper and Peierls loops evaluated at zero 
external variables are respectively given by 

Ic = - 2™W|H E E + 9c)G°_(-k) 



k>0 u„ 



tanh[e(fc)/2T(^)] 



fe>o 



e(k) 



= - -tanh[e(f)/2T] di 



(20) 



at qc = 0, where e(£) = v(£) sin^fc^ and 
I P {£) = ~ 2ttv(£)^ E E G +Ck + qp)G°- (k) 



= -Ic 



(21) 



at qp = (2feF,0). It is worth stressing that neglecting 
the dependence of Ip,c on external variables does not 
generate new momentum dependent interactions whose 
number is kept fixed along the RG flow. 

The flow equations for the local interactions then be- 
come 



digi = 

dig 2 = 
deh = 



fg9l + fl[-9i ~ 9i(92 + 92) + 9393 ], (22) 



/ S 92 + 7j/l [(93 
/ S 93 + fl[ (.92 + 92)(2ff3 + .93) 



+ h?-{gi + h)% (23) 



dm = - fgCji-, 



- (91+91X93 -| 3 )], (24) 

(25) 



where /1 = \ tanh[e(£) / 2T] . These equations differs from 
the usual scaling equations of the 1D-EG model in two 
respects. First, the rescaling for a tight-binding spectrum 
and velocity renormalization introduce linear terms; sec- 
ond, there are additional corrections coming to the cou- 
pling to momentum dependent interactions. These latter 
corrections are by far the most likely to influence the flow 
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of local couplings if not the nature of the ground state as 
we will see. 

As for the non local irrelevant interactions, the cor- 
rections due to loop contractions are small and will be 
neglected in weak coupling. From the rescaling transfor- 
mation (|13|) and the normalization of the couplings by 

nv(£), we get 

degi = — (l + Skg esc Ski + dg In v(tj) cji (26) 

for i = 1, . . .4. In the zero temperature limit, the solution 
of Eqs. ([T5|) and (f2l)|) yields the following expressions 

v(£) =v(l-— ln[2 cos 2 (^/2)]) , (27) 



and 



for the Fermi velocity and 



9i(Z) = 9i-7^e £ tan(<5^/2), 



(28) 



for the non-local couplings. The Fermi velocity is thus 
renormalized downward due to the presence of the V 
term; it reaches the value v* = v(l — ^ In 2) in the limit 
of large I. 



D. Response Functions 

To determine the nature of long-range correlations in 
the ground state, we consider the most singular response 
functions or susceptibilities, which are denoted Xu- These 
latter are obtained by adding to the £ = action an 
additional term Sur^, which consists of source fields 
linearly coupled to the composite fields O*, 



S h [r,i>] =^[M?)V^(5) + c-c], 



(29) 



Oss(q) =VjJ2 a ^-k + «)-,-aV'+,a(*), 



(32) 



k 

in the Cooper channel. Here cr^ =x ^ y ^ z (ao) are the Pauli 
(identity) matrices. 

The renormalization group transformation at the 
one-loop level, will modify Sh according to 



S h [0\0] ( 



Sh[0* , 0]i + (ShSi t 2)o,c + 



O(*)->.s O'<*) 



(34) 

where the pair fields, having zero canonical dimension, 
remain unchanged under rescaling. The last term is a 
constant oc d£zf,h*h,, that adds at each iteration and 

M M l 1 

yields the expression of the susceptibility 



™X*{fi) = 2 J ^t) z l tanh[e(l)/2T]^, (35) 

which is defined positive and evaluated in the static limit 
at q^ = (2kF,0) and (0,0) for the Peierls and Cooper 
channels, respectively. From the one-loop outer shell cor- 
rections to the linear coupling, which read 

(S h Sifl) GtC = ^t a i^[e(£)/2T]d£j2[ h AQ)9^0;(q)+c.c], 



one gets the one-loop equation for the pair vertex part 
Zp, at g°, 



where is a pair vertex renormalization factor (z M = 1 
at £ = 0). In what follows we shall examine the site 
spin density- wave (fx = SDW), bond spin density- wave 
(fi = BSDW), site charge density- wave (/i = CDW) and 
the bond order- wave (/i = BOW) susceptibilities of the 
Peierls channel; the singlet (SS) and triplet (TS) super- 
conducting susceptibilities of the Cooper channel. These 
are defined with the aid of the following expressions for 
the composite pair fields, 



O 



1 



SDW/BSDW 



(5)=5(o;„ ) ,(-?)±Ox,«,.(g)) (30) 



and 



C dw/bow(9) =2 iPoi-l) ± O (q)) 



(31) 



dt In ; 



^-tanh[e(£)/2T]. 



(36) 



For the density- wave type susceptibilities, the normalized 
couplings are given by the combinations 

.9cdw/bow = -2.9i + g 2 + g 2 T h ± 03, (37) 
3sdw/bsdw = 92 ± 9a + 92 ± §3- (38) 

The corresponding expressions for the superconducting 
susceptibilities are 



5ss/ts = Tffi - 92 T 9i ~ 92- 



(39) 



A positive value for g^ at £ — > oo signals a singularity in 
and then in x^ in that limit. 



in the Peierls channel, where 



III. RESULTS 
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Hk-q)*< ip + Ak), 



The solution of the flow equations for the pair ver- 
tices ((Ml) and the couplings (I22j|26j) in the T -> limit 
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FIG. 1: a) The phase diagram of the ID extended Hubbard 
model. The bold (thin) lines refer to the boundaries between 
the primary (secondary) phases indicated in bold (regular) 
characters. The dashed lines correspond to the boundaries of 
the phase diagram of the electron gas model in the continuum 
limit; b) zoom in the neighborhood of the U = 2V (dashed) 
line in the repulsive sector. 



leads to the determination of the most singular suscepti- 
bilities. These in turn serve to the determination of the 
dominant and subdominant phases of the model in the 
ground state. This is summarized in the one-loop phase 
diagram of Fig. [TJ as a function of weak U and V. The re- 
sults are compared with those obtained in the continuum 
limnJ-4 



A. Repulsive U 

We commence by looking at the first quadrant of the 
phase diagram, in the region surrounding the U = 2V > 
line. At the point A below the separatrix in Fig. [TJb, 




FIG. 2: a) Flow of the coupling constants 31,2,3 at the point 
A (1, 0.4) of the phase diagram in Fig.[T] b) The density-wave 
susceptibilities vs £; inset: the flow of the pair vertices de In 
for fj, = SDW, BOW and SDW. 



where U > 2V, the <?2 and gz couplings scale to strong 
repulsive values and become singular at a finite £ p , a sin- 
gularity at one- loop level that is indicative of a (Mott) 
gap in the charge sector compatible with the initial con- 
ditions satisfying the inequality gi — 2g 2 < §3- The repul- 
sive g\ coupling is marginally irrelevant and attributed 
to gapless spin degrees of freedom. The SDW response 
then develops a singularity similar to the one of BOW at 
large I ~ £ p , as shown by the behavior of zsdw and zbow 
in the inset of Fig. [5]-b. From the same Figure, however, 
the amplitude of the SDW susceptibility is larger, and 
SDW (BOW) is then taken as the dominant (subdomi- 
nant) phase in the ground state. These one-loop results 
indicate that in this region irrelevant non local couplings 
introduce no qualitative changes with respect to known 
results of the continuum theor y 1 ^ 17 . 

If we now move up to the point B in the phase diagram 
of Fig.QJb, close but below the U — 2V line, a qualitative 
change with respect to the results of the continuum limit 
emerges. While gi and 53 still scale to strong repulsive 
coupling, signaling the formation of a charge gap at a 
finite £ p (Fig. [3]-a) , the backscattering amplitude gi no 
longer scales toward zero, but extends across the g\ = 
line to then level off at a small non universal negative 
value (inset of Fig. [3]-a). According to the expressions 
in ([37]) . this change of sign of gi yields 3bow > 5sdWj 
indicating that the strongest singularity now occurs for 
the BOW response (inset of Fig. |3}b). The BOW phase 
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FIG. 3: a) Flow of the coupling constants pi, 2, 3 at point B : 
(1, 0.495) of the phase diagram in Fig.[T] b) The density- wave 
susceptibilities vs £; inset: the difference between the BOW 
and SDW flows of the pair vertices showing the dominance of 
the BOW phase. 



then becomes the dominant phase, whereas SDW closely 
follows as the secondary phase. The change of sign of </i 
takes its origin in the presence of non local couplings in 
the flow equations (I22j|24[) . Although irrelevant, these in- 
teractions push the renormalization of g\ (53) downward 
(upward) through their coupling to local variables. 

The dominance of the BOW phase becomes more pro- 
nounced as one moves up across the line U — 2V (point 
C of Fig. Hl-b). In this region, the initial local couplings 
gi and 173 are negative, but the latter interaction is still 
pushed to strong repulsive sector by non-local couplings 
(Fig. |4]-a). The BOW susceptibility then develops the 
strongest singularity with the largest amplitude (Fig. [3J- 
b). These features of the flow and the predominance of 
BOW order keep on up the BOW-CDW boundary pass- 
ing just below point D in Fig. [T}b. At that point, strong 
attractive coupling in gi and gi is occurring while g^ 
remains small and attractive (Fig.[S]-a), implying the for- 
mation of a gap in the spin sector instead of the charge. 
In these conditions, we have gcuw > 5bow, which marks 
the onset of a dominant CDW phase. The BOW order is 
subdominant and SDW correlations are non longer sin- 
gular and are strongly reduced by the presence of a spin 
gap. It is worth noting that the emergence of a spin gap 
regime on the BOW-CDW frontier is compatible with the 
results of quantum Monte Carlo simulations 7 , which find 
the onset of a Luther-Emery liquid with a spin gap at 



FIG. 4: a) Flow of the 31,2,3 couplings at C (1, 0.55) in the 
phase diagram of Fig. [T] b) The susceptibilities vs £; 



the boundary. The same analysis carried out as a func- 
tion of U allows for the delimitation of a small but finite 
fan-shape region of the weak coupling phase diagram of 
Fig. [T] where the BOW order intervenes as the ground 
state around the U — 2V line. This well known result of 
numerical calculations^— and functional RG 12 contrasts 
with the direct SDW to CDW transition predicted for 
the ID-EC modeli. 

We proceed on the analysis of the repulsive U region 
by looking at the point E, that is above the intermediate 
BOW region. In this domain, gi and 33 scale to strong 
repulsive and attractive couplings, respectively, while g\ 
is non universal and weakly attractive (Fig. [BJa), con- 
trary to what is found for the electron gas model^'i 7 -. 
The CDW singularity is stronger and accompanied by a 
weaker singularity in the BSDW response (Fig.[6}b). The 
BSDW replaces BOW as the subdominant phase over a 
finite domain of the phase diagram at V > (Fig.[l]-a). 

We turn to the point F located in the V < region 
below, but near the U — ~2V line, where qualitative 
changes with respect to the continuum limit are also 
found. In the framework of the I D-EG model, the re- 
gion below the U — —2V line is characterized by the 
conditions g\ > and g\ — 2gi > 173, respectively for 
gapless spin and charge degrees of freedom with domi- 
nant TS and subdominant SS phases^. In the presence 
of non local couplings, however, while g\ is marginally 
irrelevant, both gi and gs scale to strong repulsive cou- 
pling signaling that the charge degrees of freedom are 
still gapped (Fig. [TJ-a). Therefore the SDW phase re- 
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FIG. 5: a) Flow of the 31,2,3 couplings at D (1, 0.609) in the 
phase diagram of Fig. Q] b) The susceptibilities vs £; inset: 
the flow of the pair vertices dt lnz M for fi = CDW, BOW, and 
SDW 
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FIG. 6: a) Flow of the 31,2,3 couplings at E (1, 0.65) in the 
phase diagram of Fig. [1] b) The susceptibilities vs £. 
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FIG. 7: a) Flow of the 31,2,3 couplings at F (1, -0.55) in the 
phase diagram of Fig. [1] b) The susceptibilities vs £. 



mains dominant contrary to the 1D-EG prediction of a 
gapless TS phasei^^ (Fig- 0-b) ; the SDW incursion be- 
low the U = — 2V line expands in size as U increases as 
shown in Fig.QJa. It is worth mentioning that the result- 
ing inward bending of the TS-SDW boundary line which 
becomes more pronounced with increasing U is consistent 
with the numerical results of Nakamura 6 . 

Finally, as one moves sufficiently downward along the 
V axis, one reaches a region where gi and 33 behave 
the way marginally irrelevant variables do (Fig |8}a) , as 
shown for instance at the point G of the phase diagram of 
Fig-DJa. One then essentially recovers the behavior of the 
1D-EG model with a dominant (subdominant) power law 
singularity Xts(SS) « exp(7 TS(ss) £) for TS (SS) response 
at large I (Fig Ob) with 7ts ^ 7ss > 0. 



B. Attractive U 

We now consider the region of negative U near the 
U = 2V line. In this region, we encounter an alteration 
of the 1D-EG phase diagram boundary that is similar to 
the one discussed in the last paragraph at U = —2V > 0. 
At H in Fig. [TJb, a portion of the phase diagram with 
dominant (subdominant) TS (SS) gapless phase is lost, 
this time to the benefit of a SS phase with a spin gap. 
Strong attractive coupling in the spin sector is induced 
by non local couplings that push downward the renor- 
malization of gi (Fig. [5]-a). As for Umklapp scattering, 
it stays weakly attractive indicating that the charge sec- 
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FIG. 9: a) Flow of the 31,2,3 couplings at H (-1, -0.46) in the 
phase diagram of Fig. [T] b) The susceptibilities (logarithmic 
scale) vs £ . 



tor is gapless. The SS-TS boundary is then distorted 
inward compared to the straight line 1D-EG prediction, 
which is in fair agreement with the results of exact diag- 
onalisation by Nakamura£ The SS phase expands from 
the bent boundary with the TS phase up to the V = 
symmetry line for the transition to CDW (Fig. [TJa). We 
exemplify the SS region by the point H of the phase di- 
agram (Fig. 03-a), where the gi and 52 scale to strong 
attractive coupling for the formation of a spin gap at £ a 
(Fig. [9]-a). The SS response is the only singular response 
of the system and the whole region has no subdominant 
phase (Fig.db). 

We end the tour of the phase diagram with the second 
quadrant above the V — SS-CDW frontier at the point 
I. There, the rapid flow to strong attractive coupling for 
gi marks the onset of a spin gap at relatively small £ a 
(Fig. [TOl-a). The strong attraction for gi prevails over 
the Umklapp term, though also marginally relevant. The 
singularity of the CDW response is thus by far prevalent, 
being followed by a much weaker BOW susceptibility, 
whose subdominance is less guaranteed since it occurs in 
the strong coupling domain where the perturbative RG 
becomes less reliable. 



IV. CONCLUSION 

In this work we have proposed a generalization of the 
momentum shell renormalization group transformation 
that is applicable to ID lattice models of interacting 
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electrons. The approach has been put to the test for 
the determination of the phase diagram of the extended 
Hubbard model in weak coupling. The method discloses 
the influence of a finite number of dangerous irrelevant 
couplings on the scaling of marginal interaction terms of 
the model. Modification of scaling gives rise in some re- 
gions of the phase diagram to unexpected phases from the 
standpoint of the theory in the continuum limit. Among 
the results obtained, let us mention the incursion of BOW 
order in a finite portion of the repulsive U ~ 2V sector 
of the phase diagram, which agrees with previous results 
of numerical and functional RG methods. The approach 
is also able to capture the deformation of boundaries be- 
tween Luttinger liquid and gapped phases in the phase 
diagram of the model as found previously by exact diag- 



onalisation. 

These findings are encouraging for applications to 
other weak coupling ID or quasi- ID interacting electron 
models in which lattice details can play an important role 
in the properties of correlations at long distance. 
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